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Abstract. We use publicly available l/-band imaging data from the wide-angle surveys conducted by the ESO 
Imaging Survey project (EIS) to further investigate the nature of the EIS galaxy cluster candidates. These were 
originally identified by applying a matched-filter algorithm which used positional and photometric data of the 
galaxy sample extracted from the J-band survey images. In this paper, we apply the same technique to the galaxy 
sample extracted from V-band data and compare the new cluster detections with the original ones. We find that 
~ 75% of the low-redshift cluster candidates {z < 0.6) are detected in both passbands and their estimated redshifts 
show good agreement with the scatter in the redshift differences being consistent with the estimated errors of the 
method. For the "robust" 7-band detections the matching frequency approaches ~ 85%. We also use the available 
{V — I) color to search for the red sequence of early-type galaxies observed in rich clusters over a broad range 
of redshifts. This is done by searching for a simultaneous overdensity in the three-dimensional color-projected 
distance space. We find significant overdensities for ~ 75% of the "robust" candidates with zi < 0.6. We find 
good agreement between the characteristic color associated to the detected "red sequence" and that predicted 
by passive evolution galaxy models for ellipticals at the redshift estimated by the matched-filter. The results 
presented in this paper show the usefulness of color data, even of two-band data, to both tentatively confirm 
cluster candidates and to select possible cluster members for spectroscopic observations. Based on the present 
results, we estimate that ~ 150 EIS clusters with zi < 0.6 are real, making it one of the largest samples of galaxy 
clusters in this redshift range currently available in the southern hemisphere. 
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1. Introduction 

Clusters of galaxies are prime targets in observational cos- 
mology for several reasons among which: they are the 
largest gravitationally bound structures observed in the 
universe and as such are useful tracers of the large-scale 



structure; Llielr abuiidauw; as a iuiicLiuu oi redsmiL can 
strong cousLraluLs on Llie cusmulugical parameters; sUid- 



rently not available, with most of the existing cluster cat- 
alogs being either limited in number or redshift coverage. 
In recent years, several attempts have been made to con- 
duct systematic searches for more distant systems both 
in optical/infrared and X-ray utilizing a large variety of 
detection techniques (e.g. Gunn et al. 1986; Lidman & 
Peterson 19961; iLumsden et al. 19921 [Postman et al. 19961; 



^aritsky et al. 1997 ; Gonzalez et al. 2001; pladders fc 



ies of the pruperLles uf cluster galaxies can pruviue msi] 
on theJi sLai-ruimaLiuu and meigiug hisLuiy and set valu- 
able constraints to galaxy formation and evolution mod- 
els. However, such studies are only possible using statis- 
tically representative samples of clusters spanning a large 
range of redshifts. Unfortunately, such samples are cur- 



piovide maifehl yee 2000|). Other jTomising methods include the use of 
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mass reconstruction maps built from weak lensing signals 
and the Sunyaev-Zeld'ovich effect to conduct blind cluster 
searches over wide areas. 

While the number of candidates has significantly in- 
creased, especially from optical/infrared searches, the 
number of confirmed and well-studied cases with z ^ 



1 remains relatively small (e.g. Postman et al. 2001) 
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Therefore, the present samples fall far short of what is 
desirable for: exploring the relationship between systems 
selected using different detection algorithms, a more com- 
plete characterization of cluster properties and a better 
understanding of the evolution of clusters and their galaxy 
population. The major obstacles for enlarging the sam- 
ple of well-studied cluster galaxies have been a combina- 
tion of several factors which include the reliability of clus- 
ter candidates, interloper contamination and instrumen- 
tal limitations. The reliability issue has in general favored 
follow-ups of X-ray selected samples, which are currently 
small and unlikely to grow significantly in the near fu- 
ture. Therefore, to immediately benefit from the available 
large aperture telescopes and multi-object spectrographs 
for galaxy cluster research, for the time being one must 
rely on optical/infrared surveys such as those being carried 
out by the ESO Imaging Su rvey (EIS) project. Using the 



combined with the results of ongoing spectroscopic obser- 
vations being conducted at the 3.6m telescope at La Silla 
(e.g. Ramella et al. 200C) and at the VLT of color-selected 



original /-band survey data (Nonino et al. 199E; Prandoni 
et al. 1 



j99| ; Benoist et al. 199£ ) a catalo g of cluster candi- 
dates was compiled ( Olsen et al. 1999a| ,|b|; ^codeggio et al. 



1 999|)~^omprising > 300 candidates in the redshift range 
0.2 < z < 1.3 within an area of ~ 15 square degrees, 
until recently the largest publicly available sample of its 
kind in the southern hemisphere. However, candidate clus- 
ters based on a single passband are normally viewed with 
skepticism because of the known problems caused by pro- 
jection effects. In order to address this concern, in the 
present paper we use the publicly available V^-band data 
to demonstrate that by using color information one can 
provide independent supporting evidence of the reality of 
a large fraction of the candidate clusters. 

We use the available V data in two complem entary 



ways. First, w e extend the analysis carried out by Olsen 
et al. 1999b| ) and treat the data from the V^-survey as 



an independent data set to which we apply the matched- 
filter algorithm and test the reproducibility of the detec- 
tions. The F-detections are then cross-compared to those 
obtained from the analysis of the /-band data. Second, 
we use the {V — I) color to search for a simultaneous 
spatial and color concentration at the location of the / 
detections. Such concentrations are interpreted as being 
associated with the red sequence observed in rich clusters 
over a broa d range of redshifts ( e.g. Aragon-Salamanca 

[Lubin 1996| ). This red 



et al. :_993: Stanford et al. 199 



sequence is interpreted as being an indication that rich 
clusters contain a core population of passively evolving el- 
liptical galaxies which are old and coeval. This method is 
similar to that described by [Gladders fc Yee (2000 ) in the 



so-called red sequence cluster survey. 

The y-band data restricts the analysis primarily to 
cluster candidates with z < 0.6 since at larger redshifts 
ellipticals dropout from the sample as the 4000A break 
moves past the y-band. However, the analysis presented 
here illustrates the methodology that will be employed in 
forthcoming papers to investigate high-redshift candidates 
using the multi-color optical/infrared {R and JKs) imag- 
ing data assembled for a large number of EIS candidates 
with z > 0.5 (~ 110 and ~ 40). These data will also be 



systems in different redshift intervals. Preliminary results 
based on an admittedly small sample suggest that a signif- 
icant fraction of these color-selected systems are confirmed 
spectroscopically. 

This paper is organized as follows. Sect. ^ briefly de- 
scribes the data used in the present analysis. Sect. |^ 
presents the results obtained from the comparison of 
matched-filter detections in V and /, from which a first 
estimate of the confirmation rate is derived. Sect. ^ de- 
scribes the methodology employed to detect early-type 
cluster members and to measure the typical color asso- 
ciated to these objects. This information is used to obtain 
independent estimates of the expected confirmation rate 
and the redshifts assigned to them. A discussion of the 
results is presented in Sect. ^. A brief summary of our 
results is given in Sect. I 

2. The Data 



In the present analysis we use the catalog of EIS can- 
didate clusters compiled by Olsen et al. (19994 [b|) and 
IScodeggio et al. (1999D comprising > 300 candidates over 
four patches of the sky covered by the EIS- WIDE /-band 
survey conducted at the NTT. The candidates were iden- 
tified from the galaxy distribution extracted from these 
images and presented in Nonino et al. (1999), Prandon: 
|et al. (1999| ) and [Benoist et al. (1999| ). 

In this paper the /-data are complemented with those 
taken in V. The F-data consists of two sets of images. 
One taken at the NTT with EMMI, as part of the orig- 
inal EIS- WIDE survey which covered two square degrees 
(Patches A and B). The other consists of images taken 
at the MPG/ESO 2.2 m telescope with ESO's Wide-Field 
Imager (WFI) as part of the so-called Pilot Survey. The 
latter set covers a total area of 12 square degrees in V- 
band, out of which 6 square degrees are publicly available. 
Galaxy catalogs have been extracted from these ^-images 
and associated with those in / to produce {V—I) color cat- 
alogs utilized below, covering a total area of 8 square de- 
grees, or about half the area of the /-survey. Within this 
area there are 168 /-band cluster candidates, correspond- 
ing to 55% of the entire EIS candidate sample. 

Fig. H compares the galaxy counts obtained from the 
V^-band catalogs extracted from the images taken using 
different telescope/detector setups. Even though reaching 
different limiting magnitudes, the counts are consistent 
with each other and in good agreement with those of other 
authors (Arnouts et al. 1995). This demonstrates the uni- 
formity of the galaxy catalogs, which are not distinguished 
in the analysis below. The 5g-limiti ng magnitudes are es- 
timated to be ^ 24 and V - 24.6 (|da Costa et al. 2000|) 
for the EIS- WIDE (patches A and B) and Pilot Survey 
data (patch D), respectively. Therefore, in the analysis 
presented below we only consider galaxies brighter than 
V = 24. 
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Fig. 1. Galaxy number counts derived from the V^-band 
images obtained for EIS Patches A (dashed-dotted hne) 
and B (dashed line) and the Pilot Survey (Patch D, dotted 
line). For comparison the counts of Arnouts et al. (1995) 
(open circles) are also shown. 



3. Searching for clusters in l/-band 

3.1. Matched-filter Detections 

The EIS cluster search was using the matched-filter al- 



gorithm originally proposed by Postman ct al. (1996) 
and used to build the Palomar Distant Cluster Survey 
(PDCS). This method uses a filter built from the luminos- 
ity function and radial profile of nearby clusters. At each 
position of the survey area the probability of the presence 
of a cluster matching the filter is computed. By fine-tuning 
the filter to match the apparent magnitudes and extent of 
the cluster at a series of redshifts (0.2 < z < 1.3) a like- 
lihood curve is constructed. The cluster candidate is then 
assigned the redshift for which the maximum likelihood is 
found. 

In contrast to the PDCS the EIS cluster catalog was 
built exclusively from matched-filter detections in the I- 
band for which two sets of images, contiguously cover- 
ing the regions of interest, were available. These sets were 
treated separately to construct two galaxy catalogs which 
were called "odd" and "even" catalogs, to indicate the im- 
age mosaic from which they were extracted. The matched- 
filter algorithm was run on each catalog and the result- 
ing cluster detections were compared. The detections were 
considered to be: "robust" if they were detected at 4ct in 
one or both catalogs, or if they were detected at the 3a- 
level in both the "even" and "odd" catalogs, and "poor" if 
they were detected at the Scr-level in only one of the cat- 



alogs. These definitions were introduced by Olsen et al 



1999a) (see also Scodcggio et al. 1999) and are used be- 



low. 

The V images available represent another sampling of 
the galaxy distribution in the region and as such can be 
used to define another set of cluster candidates. In apply- 
ing the matched-filter procedure to the galaxy catalogs 
extracted from the F-band images, we use the same clus- 
ter model as in /. This model assumes a Schechter lumi- 



nosity function and a radial distribution of galaxies given 
by a truncated Hubble profile (for more details see Olsen 
et al. 1999a). Following Postman ct al. (1996) we adopt 
the Schechter parameters My = —21.02 and a faint end 
slope of a — —1.1, valid for local clusters. For the radial 
profile we adopt a core radius of Vc — 100/i^^kpc and a 
cut-off radius of Vco = lft.~^Mpc. We use a cosmological 
model with Hq = 75km/s/Mpc, fig = 1.0. 

In the discussion presented below it is important to 
have a rough idea of the redshift limit for cluster de- 
tection in each of the passbands. This can be obtained 
from the limiting magnitudes of the input galaxy cata- 
logs (/lim ^ 23 and ~ 24, Vega system), by assuming 
that the cluster population is formed predominantly by 
early-type galaxies, and by assuming a particular model 
for the evolution of the stellar population. The maximum 
contrast of a cluster over the field population occurs in 
the range of magnitudes m* ± Imag, where m* is the 
apparent Schechter magnitude of the cluster. Therefore, 
the limiting redshift corresponds approximately to that 
for which m* -I- 1 is comparable to the limiting magnitude 
of the galaxy catalog. Assuming a no-evolution model this 
would, for the data considered here, imply a redshift limit 
of z 0.6 and z ~ 1.1 for the V- and /-candidates, respec- 
tively. If one instead assumes a passive-evolution model 
the redshift limits would become slightly larger. As it will 
be seen below, z ^ 0.6 as determined for the y-band is 
indeed the practical limit of the present analysis. 



3.2. Comparison between V- and I-detections 

The cluster candidates identified by applying the 
matched-filter algorithm to the galaxy catalogs extracted 
from the V images are compared to the /-band candidates. 
Below we only consider those which match /-candidates, 
even though there are also ^-candidates without corre- 
sponding /-detections (^ 40, which will be discussed else- 
where) . Detections in the different passbands arc matched 
by considering only positional coincidence. For this we re- 
quire the angular separation to correspond to a projected 
separation of less than 0.5/i~^Mpc at the redshift esti- 
mated for the /-candidate (z/), regardless of the matched- 
filter redshift estimate, zy, based on the V catalogs. In 
the following we investigate the frequency of matches and 
compare the V- and /-estimated redshifts. 

Fig. 1^ summarizes the results of the cross-comparison 
between V- and /-candidates. The top panels show the 
redshift distribution of the matched candidates compared 
with the redshift distribution of the sub-sample of EIS /- 
candidates located within the area covered by V. The left 
panels refer to "robust" candidates (90 candidates) and 
the right to the "poor" ones (78 candidates). The overall 
matching rate is 61% (102 out of 168 /-band candidates), 
and varies with redshift as shown in the lower panels of 
Fig. |. 

From Fig. ||one sees that for z/ < 0.6 most clusters are 
detected in both passbands. In this redshift interval the 
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Fig. 2. The matching rate as function of redshift for the "robust" (left cohinin) and "poor" (right column) candidates 
as explained in the text. The upper panel shows the redshift distribution for all the candidates (dashed histogram) and 
those with counterparts in the F-band (filled histogram). The lower panel shows the respective matching frequency 
as function of redshift which is interpreted as an estimate of the confirmation rate of the original /-detections. 



matching rate (hereafter interpreted as an estimate of the 
confirmation rate) for the whole sample is ~ 75% (88 out 
of 118). Additionally restricting the sample to "robust" 
candidates, the confirmation rate is higher, reaching ~ 
82% (69 out of 81). It is also interesting to note that about 
50% of the "poor" candidates in this redshift range arc 
associated with ^-detections, showing an above average 
rate at zi = 0.3. This could be due to low-redshift poor 
clusters detected at a low significance level in a single- 
passband. 

At larger redshifts {zj > 0.7) only 14 (~ 26%) of the 54 
7-band cluster candidates are matched by a ^-detection. 
This lower confirmation rate is not surprising considering 
that for systems with z > 0.6 the 4000A break lies red- 
wards of the F-band. Note, however, that we find match- 
ing y-detections over a broad range of redshifts including 
some very high-redshift systems {zj > 1), especially in the 
"poor" sample. In order to investigate the nature of these 
matches the 14 cases of clusters with zj > 0.6 and with 
matched V detections were visually inspected. In five cases 



{zi = 0.7,0.7,0.9, 1.1, 1.3) the F-detection is clearly asso- 
ciated with a foreground concentration unrelated to the 
/-band cluster candidate, even though the /-band image 
does appear to include a cluster candidate. In one case 
{zi = 1.3) the F-image contains a satellite track which 
leads to spurious galaxy detections forming an artificial 
concentration of objects in the V catalog. For the remain- 
ing eight cases the ^-band images do show concentrations 
of galaxies consistent with those seen in /. These eight 
candidates - four robust and four poor /-detections - have 
redshifts in the range 0.7-1.1. As mentioned above, while 
early-type galaxies arc not expected to be detected in this 
redshift range, the cluster detection based on the matched 
filter technique takes into account the entire galaxy pop- 
ulation. Therefore, this makes this method less sensitive 
than others (see below) to the fading in the F-band of 
the early- type galaxies. One should also keep in mind the 
large uncertainties in the redshift estimates at high red- 
shift, which could lead to the matching y/-detections with 
large differences in redshift estimates. However, as shown 
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the cluster seems real in both bands. However, the orig- 
inal /-detection had two different redshifts in the "even" 
{zj ~ 0.4) and "odd" (z/ — 0.2) catalogs. The latter was 
included in the EIS cluster candidate catalog, because it 
had a higher significance. Clearly, for this cluster zj is 
uncertain, but there is little doubt about the reality of 
the cluster. In the second case the redshift discrepancy is 
probably due to the superposition of a low redshift sys- 
tem detected in V and a more distant cluster detected 
in /. Finally, we emphasize that the above results also 
demonstrate the reliability of the redshift estimates and 
can be used for selecting sub-samples in different redshift 
intervals, if an error of about 0.15 is properly taken into 
account. 



Fig. 3. The distribution of differences of redshifts esti- 
mated from the V and / matched-filter identifications for 
the entire sample (dashed line) and the "robust" candi- 
dates (solid line). 



below this does not seem to be the case. Of course, there 
is still the possibility of the superpositions of clusters at 
different redshifts with / detecting a distant richer sys- 
tem, while V detecting a foreground one. Deciding among 
these possibilities will have to await the results of spec- 
troscopic observations. Note that a comparable fraction 
of VI detecti ons in the redshift in terval bei ng considere d 
was found by Postman et al. (199(S| ) (see also Lubin 199^ ) . 



3.3. Redshift Comparison 

Besides confirming /-candidates, the V-detections can 
also be used to assess the reliability of the original clus- 
ter redshift estimates based on the application of the 
matched-filter algorithm to the /-band galaxy catalog. 
Fig. ^ shows the distribution of redshift differences es- 
timated by the matched-filter algorithm for associated V- 
and /-detections. The figure shows this distribution for 
the whole sample, considering matches over the entire red- 
shift range, but eliminating the six spurious matches dis- 
cussed in the previous section. Also shown is the distri- 
bution of redshift differences for the sample of "robust" 
candidates. Taking the sample as a whole, we find a mean 
difference of (zy — zj) — —0.03 and a scatter of 0.12. 
For the "robust" candidates the mean difference is sig- 
nificantly smaller (zy — zj) — —0.003 with a scatter of 
0.11. Interestingly, all the eight cases with matched VI 
detections and redshifts zi >0.6, discussed in the pre- 
vious section, have consistent redshifts within the esti- 
mated errors, giving further credence to their confirma- 
tion. Note that two of the "robust" cases are outliers - 
one with (zy — zj) = 0.5, zy = 0.7 and zj = 0.2; an- 
other with {zy — zi) — —0.4, Zy — 0.2 and zj — 0.6. To 
investigate these cases we visually inspected the i mage s 
and the "even" and "odd" detection fists (see Sect. 3.1). 
From this inspection we conclude that: in the first case 



4. Searching for cluster early-types 

Confirmation of candidate clusters by matching V and / 
relies only on the fact that the V images and the derived 
galaxy catalogs may suffer different effects than those in 
/-band. Thus, the reproducibility of the results obtained 
from data in different passbands serves to corroborate the 
original findings, thereby increasing the significance of the 
detections. This test, however, relies on the same underly- 
ing assumptions of the matched-filter method - clusters are 
nearly spherical, have the same projected density profile 
and are characterized by a Schechter luminosity function, 
similar to the one observed in local clusters. 

An alternative way of searching or checking the reality 
of clusters is to use the additional information provided by 
the {V — I) color. One way is to search, at the location of 
the / detections, for early-type galaxies which are known 
to form a tight color-magnitude relation and to populate 
the central regions of galaxy clusters. Another is to con- 
duct blind searches for clustering of early-type galaxies 
with colors consistent with their being at different redshift 
intervals. These methods make no assumption about the 
global properties of clusters but have to rely on galaxy evo- 
lution models. These color-based methods are complemen- 
tary to the matched-filter approach. In addition, splitting 
galaxy catalogs in color intervals leads to an increase in 
the contrast between the cluster and the background. This 
enhanced contrast can be particularly important for the 
confirmation of poor cluster candidates. Consistency be- 
tween the two approaches, for reasonable models of galaxy 
evolution, can thus lend further support to the reality of 
a cluster candidate. Additionally, an important advantage 
of color-based methods is that they provide a sample of 
candidate early-type cluster members which can be used 
for follow-up spectroscopic observations. 

4.1. Color-slicing procedure 

In order to detect clustering of early-type galaxies associ- 
ated to the /-detections we search for simultaneous over- 
densities in color and projected density distribution. We 
select galaxies within a region of 3 x 3ft-~^Mpc^ and within 
a magnitude interval of m* ± 2, where m* is the expected 
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Fig. 4. Six density maps for one cluster candidate. The upper and lower row each correspond to the first three maps 
of each of the subsets of seven density maps. The remaining eight slices do not contain relevant information in this 
particular case. The two sets are shifted by half the width of a color slice to make sure that a cluster is not split 
between two slices, whereby it would be less likely to be detected. The contours mark the regions with densities in 
excess of the mean density and are separated by one a. The shaded areas mark the significant pixels. The scale on the 
axes is arcmin. In each panel the adopted pixelscale is marked as the filled square in the lower right corner. 



apparent Schechter magnitude, around the /-candidate 
position and split the galaxy sample in color intervals 
A(y - /)=0.3, covering the range 1.1 < (V" - /) < 3.35 
corresponding to the redshift range < z < 1.2. For each 
color interval we construct a smoot h density map using 
an adaptive kernel smoothing (see Merritt fc Tremblay 
1994[)~| The pixel scale of the density map scales with 
the estimated redshift of the candidate considered. In the 
present analysis it is chosen to be 75/i~^kpc, which roughly 
corresponds to the cluster model core radius. We select 
significant pixels at the 3a confidence level relative to the 
field galaxy density distribution in the color interval be- 
ing considered. The density peaks are then identified in 
the spatial distribution of these pixels. The background is 
computed averaging the smoothed density distributions of 
60 randomly selected areas chosen within the EIS patches. 
If a significant peak lies within a radius of 375/i"^kpc (cor- 
responding to 5 pixels) from the nominal cluster position 
we consider it to be a detection of the cluster candidate 
red sequence. We interpret these detections as a confirma- 
tion of the /-detection. The search radius adopted allows 
for possible offsets of early-type galaxies relative to the 
cluster center as defined by the matched-filter position. 
In cases where significant pixels are found in more than 
one color interval the color assigned to the cluster red se- 
quence is that of the slice in which the largest detected 



overdensity is found. In the case of superpositions in the 
cluster fields we manually assign the color that is closest 
to that predicted by galaxy evolution models. 

Fig. U illustrates the results of the color-slicing analy- 
sis showing two sets of density maps obtained in different 
color intervals. In this particular case, significant pixels 
are detected in two different slices. The highest density is 
obtained in the slice of 1.10 < {V — I) < 1.40, thus a color 
{V — 1)=1.25 is assigned to the cluster. More details about 
the color slicing procedure adopted here can be found in 
plsen (200C ). Here it suffices to say that most of the orig- 
inal /-candidates with zj < 0.4 exhibit well-defined red- 
sequences which are clearly visible in the (^ — /) x / color- 
magnitude diagram. Furthermore, the characteristic color 
of the observed red sequence is, in general, consistent with 
that derived by the automatic procedure described above. 
For clusters with zj > 0.4 the data are not sufficiently 
deep for the red-sequence to be clearly visible for individ- 
ual clusters {e.g., Lubin 1996) and we must rely on the 
color-slicing detections. 

4.2. Detection of the red sequence 

Applying the color-slicing technique described above to 
the /-band candidates we find significant detections for 94 
out of 168 56%) candidates, a value comparable to that 
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Fig. 5. Each panel shows the histogram of {V — /)-colors assigned to the cluster candidates with estimated redshift as 
indicated above each panel. The dashed bars indicate the {V — /)-interval corresponding to a redshift range of z ± 0.1 
corresponding to the estimated uncertainty. The solid bars mark the {V ~ /)-interval centered at z — 0.1 as discussed 
in the text. The colors are estimated by a passive evolution model for the stellar population as described in the text. 



previously found for matched V- and /-detections. The 
redshifts of the candidates cover the range from zj — 0.2 
to Zi = 1.1. Of the candidates with significant detections 
69 were also detected by the matched filter in V, which 
correspond to 68% of the 102 candidates detected in both 
V and / by the matched filter. When considering only the 
cases with zj < 0.6 for which the y-band filter is most 
suitable, we find that 79% of the V and / matched-filter 
detections are confirmed by the color slicing. In order to 
investigate the colors assigned to each candidate we show 
in Fig. H the histograms of the colors assigned to candi- 
dates in each redshift bin. Redshift bins without color-slice 
detections are not shown. The intervals shown on the top 
of each panel indicate the uncertainty in the (V — /)-color 
due to errors in the redshift estimates. In this plot the 
(V — /)-color interval is computed using a model for a 
passively evolving stellar population for a redshift range 



of z ± 0.1 (dashed line). The adopted model assumes a 
Salpeter initial mass function and a burst of star forma- 
tion lasting less than 1 Gyr and an epoch of formation 
corresponding to Zf = 10 (see further discussion below). 
One immediately notices that the predicted colors are red- 
der than those obtained for the cluster candidates. This is 
consistent with results based on simulated data showing 
that the matched -filter procedure overestimates the red- 
shift ( Olsen 2GGC| ). In Fig. || we also show the color interval 
(solid line) predicted by the model taking this zero-point 
offset into account (see below). Note that the {V — /)- 
uncertainty becomes smaller at z ~ 0.5 — 0.6, because at 
these redshifts the evolution of the color with redshift is 
slow. 

From the figure one finds that in the redshift range 
z/ ^ 0.2—0.6 the color assigned to the cluster red-sequence 
form a well-defined distribution consistent with model pre- 
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Fig. 6. The panels to the left show the results for the sample of "robust" candidates and the right column for the 
"poor" ones. Upper panels: The redshift distribution of clusters detected by the color slicing method (dark histogram) 
compared to the distribution for all candidates with V- and /-data (dashed line). Lower panels: The confirmation rate 
computed as the number of candidates with a color slice detection (the dark area in the upper panel) relative to the 
total number of candidates with V- and /-data (the dashed line in the upper panel.) 



dictions if, as mentioned above, one assumes that the red- 
shift estimates, based on the matched-filter, are indeed 
overestimated. However, some outliers are clearly present, 
their number increasing with redshift. 

In the redshift interval zj ^ 0.5 — 0.6. wc find seven 
blue and three red outliers. In order to investigate their 
nature the images and density maps for these detections 
were individually inspected. In all cases the position of 
the detections correspond well with that of the cluster 
candidate. Moreover, for all of the blue outliers but one, 
the color assigned by the automatic procedure is associ- 
ated to foreground detections. However, in most of these 
cases redder peaks are present but at lower significance, 
preventing them from being considered detections by the 
automatic procedure. These redder peaks are probably re- 
lated to the cluster red sequence, and thus represent the 
real color of the early-type members. In this paper, we take 
a conservative stance and do not consider these cases as 



true detections and discard them from any further analy- 
sis. Therefore, we conclude that the increase in the number 
of blue outliers is due primarily to the fading of early-type 
galaxies as the 4000A break moves through the V-band fil- 
ter and that the size of the search radius adopted has no 
effect on the number of outliers. This effect is unavoid- 
able and its onset depends both on the redshift of the 
cluster and on its nature {e.g., composition, number of 
bright early- type galaxies). There are two possible ways 
of further characterizing these clusters. One is to consider 
the spatial distribution of F-dropouts, the other to use 
/i-band data collected for this project. This will be con- 
sidered in a separate paper. It is worth noting that in the 
redshift interval zj ~ 0.5 — 0.6, we also find three red 
outliers, all of which are consistent with the interpreta- 
tion that the red sequence of the cluster has indeed been 
detected. 
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At zi > 0.7 we find five candidates exhibiting a color 
slice detection. Since at these redshifts the 4000A break 
has completely moved through the V-band filter we expect 
that most of the early-type galaxies have dropped out from 
the sample and that no red sequence can be detected. In 
fact, from the inspection of the images and density maps 
we conclude that all the detections are due to foreground 
concentrations. This conclusion is also valid for the two 
cases with matching V and / matched-filter detections 
(zy=0.7, zi = 0.8; zy=0.9, zj = 1.0). While these systems 
are detected by the matched-filter technique, which relies 
on the overdensity of the entire cluster galaxy population, 
they would not be detected based on the the color-slicing 
method for reasons mentioned above. 

In Fig. 1^ we compare the redshift distributions of the 
/-detections with those for which a red sequence has been 
detected. As before, we show the sub-samples compris- 
ing "robust" and "poor" candidates, separately. The lower 
panels show the estimated confirmation rate as a func- 
tion of the matched-filter redshift zj derived for the I- 
candidates. 

For z/ < 0.6 we find 80 clusters confirmed out of 118 
(~ 68%). Restricting the sample to only "robust" can- 
didates we confirm 60 out of 81 candidates 74%). In 
this redshift range we find that the confirmation rate is 
nearly constant, except for the outermost bins, both show- 
ing lower confirmation rates. We point out that for very 
nearby systems zj ~ 0.2, the red-sequence may not be de- 
tected automatically due to the large angular size of the 
clusters. For the zj ~ 0.6 the lower rate possibly refiects 
the fading of early-type galaxies due to the redward shift 
of the 4000A break. Finally, note that for Zj < 0.6, ~ 50% 
of the "poor" candidates are confirmed by the color-slicing 
method, in agreement with our earlier findings based on 
matched V- and /-detections. 

Clearly, the color-slicing method, which is based on 
the identification of elliptical galaxies, is more sensitive to 
their fading in the l^-band than the method that relies on 
the detection of matched V- and /-detections which de- 
pend on the galaxy population as a whole. This restricts 
the effectiveness of the identification of red sequences up 
to roughly z 0.6 and thus explains the sharp drop ob- 
served in the confirmation rate beyond z/=0.6. This effect 
combined with variations in individual cluster properties, 
the uncertainties in the redshift estimate and colors, and 
projection effects might also explain the significant scatter 
in the estimated color (discussed in the next section) of 
the red sequence for systems with zj > 0.5. 

4.3. Colors of early-type cluster members 

The color-slicing method provides an estimate of the color 
of the detected red sequence of a cluster candidate, which 
can be used to check the consistency with that predicted 
by a particular galaxy evolution model at the estimated 
redshift zj. This can be used either to constrain the epoch 
of formation of ellipticals or inversely to convert a mea- 




Fig. 7. The relation between the assigned {V — /)-color 
and estimated redshifts for the clusters that are considered 
confirmed. The filled circles mark the mean color in each 
redshift bin with Icr-error bars. The dashed line is the 
color evolution for a non-evolving elliptical galaxy (only 
applying the appropriate K-corrections for the redshifted 
galaxy spectrum), and the solid line indicate the color 
evolution for a passively evolving galaxy with Zf = 10 



(Bruzual & Chariot 1993). The dotted line shows the locus 
of the passive evolution shifted by Az = +0.1. The dot- 
dashed line shows the locus of the passive evolution of a 
galaxy formed at Zf = 2 shifted by Az = +0.1. 



sured color to an independent redshift estimate, as dis- 
cussed below. 

Fig. ^ shows the mean colors and the 1 a error bars of 
the cluster red sequence, as assigned by the color-slicing 
method, as a function of z/. For comparison, we also show 
the colors pre dicted for non-evolv ing (dashed line, using 
spectra from Kinney et al. 19'9^) and passively evolving 
(solid line, Bruzual fc Chariot 1993| ) ellipticals formed 
at Zf = 10 as a function of z. As also seen in Fig. ^ one 
immediately finds that the matched-filter estimated red- 
shifts overestimate the true redshift by about 0.1 in order 
to match the colors of present-day ellipticals. Instead of 
correcting the z/ values we shift the curve corresponding 
to the passive evolution model to take into account this 
zero-point offset (dotted line). We also include a model 
for passively evolving galaxies, adopting a formation red- 
shift of Zf = 2, already corrected for the zero-point offset 
mentioned above (dot-dashed line). It is worth pointing 
out that the model curves are neither very sensitive to the 
adopted cosmological model nor do the predicted colors of 
ellipticals vary significantly in the redshift range consid- 
ered here. Currently, several studies put a lower limit on 
the formation epoch of ellipticals of Zf ~ 2. From the figure 
it can be seen that we find good agreement between the 
measured mean and predicted colors. The slightly larger 
scatter at z > 0.4 is probably due to larger errors in the 
matched-filter estimated errors and to the effect of the 
4000A break being shifted into the spectral region of the 
V^-band filter. 
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Fig. 8. The distribution of redshift differences between 
color estimated redshifts and the matched filter estimated 
redshifts. The color estimated redshifts are based on the 
Bruzual fc Chariot (199^ ) passive evolution scenario ex- 
plained in the text. The dashed line shows the results for 
the sample as a whole while the solid line show the "ro- 
bust" sample results. 

4.4. Redshift comparison 

Given the overall agreement between model predicted col- 
ors and those estimated from the color slicing procedure, 
the (V — J)-color assigned to the red sequence of clus- 
ter early-types can, alternatively, be used to obtain an 
independent redshift estimate for those clusters that are 
confirmed. Using the passive evolution model introduced 
earlier, we convert red-sequence colors into redshifts. A 
comparison between these redshifts and those based on 
the matched-filter is shown in Fig. |[ where the distribu- 
tion of the differences between these redshift estimates is 
presented. We find that in the mean the matched filter 
overestimates the redshift by Az ~ -1-0.11 with a scatter 
of <t(Az) ~ 0.12 both for the "robust" and the entire sam- 



nosity distribution, while the color-slicing method relies on 
the fact that in real clusters early-types are located pref- 
erentially in the central regions and form a tight color- 
magnitude relation. Similarly, the redshift estimates are 
independent, one relying primarily on magnitudes while 
the other on colors. 

Our main results can be summarized as follows. For 
z < 0.6 the overall confirmation rate by either method 
is high ~ 65% for the color-slicing and ^ 79% for the 
matched-filter confirmation. Considering only the "ro- 
bust" /-band candidates this rate becomes ~ 75% for the 
color slicing and ^ 85% for the matched-filter technique. 
Moreover, about 57% of the candidates are confirmed by 
these independent methods. This rate increases to 66% 
when only "robust" detections are considered. These re- 
sults enable us to adequately rank the candidate clusters 
for follow-up observations. For candidates at higher red- 
shifts {z > 0.7), nearly all "poor" detections, it is occa- 
sionally possible to tentatively confirm their existence by 
matching 1/ /-detections (~ 26%). 

For the confirmed clusters we find good agreement be- 
tween the redshifts estimated by the different methods 
with a scatter comparable to that expected for redshifts es- 
timated by the matched-filter procedure. We do find that 
the matched-filter estimated redshifts may be biased high 
by Az 0.1. Correcting the redshifts by this amount we 
find a very good match between the mean color of early- 
type galaxies at different redshifts and those predicted by 
the passive evolution models, consistent with the findings 



of other authors (e.g. Lubin 1996; Stanford et al. 1998 



Aragon-Salamanca et al. 1993). These results show that, 



in general the matched filter redshifts can be used as a 
rough guide for a redshift selected sample. 

Finally, our results are also consistent with the blind 
search of clusters base d on the clustering o f color-selected 
galaxies carried out by Bcnoist et al. (2001 ). In contrast to 
our analysis, which considers the clustering of red galaxies 
at the position of the matched-filter candidates, Benoist 



to the red outliers m ig. |^. 



5. Discussion 

In this paper we have argued that direct multi-band imag- 
ing data can be used to provide supporting evidence for 
the reality of candidate clusters of galaxies derived from 
the application of the matched-filter algorithm to a single 
band galaxy catalog. To illustrate this point we have, in 
this paper, used publicly available VI data to compare 
V and / detections and to use the {V — I) color to iden- 
tify simultaneous concentrations in space and color, pos- 
sibly reflecting the existence of galaxies along the color- 
magnitude relation observed in real clusters. The nature 
of the two methods is significantly different both regarding 
their underlying assumptions, the nature and methodol- 
ogy of detections and the resulting estimates of the red- 
shifts. The matched filter assumes that clusters are spher- 
ical concentrations following a specific density and lumi- 



Ple- Tlje outhcrs with (ZsUgc - zi) ^ 0.3 - 0.4 corresponds et al. (200l|) builds a sample of cluster candidates from the 



clustering of galaxies split according to the expected colors 
of early-type galaxies at different redshifts. The analysis 
considers wide-angular regions and the algorithm used in 
identifying significant peaks is also different. Nevertheless, 
there is a good agreement between the cluster sample de- 
rived and those identified by the matched-filter. 

6. Summary 

This paper presents the first results of a systematic follow- 
up program using direct imaging data to confirm EIS 
galaxy clusters. We have used publicly available V- and 
/-band data over ~ 8 square degrees to investigate the 
reliability of the original cluster candidates. To this end 
we have used two largely independent methods based on: 
a direct comparison between V- and /-matched-filter de- 
tections; second, a search for early-type cluster members 
at the position of the /-detections using the available 
{V — I) color galaxy catalogs. In the redshift interval 
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z < 0.6, for which the V-band filter is most suitable, both 
methods provide supporting evidence that most clusters 
(~ 95 systems or ~ 80% of the sample), originally identi- 
fied from the analysis of the /-band galaxy catalogs, are 
real. Interpreting these results as a true estimate of the 
expected yield of real clusters, we speculate; that when all 
the available l/-band data are fully reduced and analyzed 
a sample comprising > 150 clusters with ^; < 0.6 will be 
available in the southern hemisphere for statistical stud- 
ies. We also find that in this redshift range the estimated 
redshifts seem to be robust within Az ~ 0.1 and can be 
used to select clusters in different redshift bins. To extend 
the present analysis to higher redshifts it is necessary to 
use redder passbands. For this purpose we have assembled 
R- and Jii's-band data for ~ 110 and ~ 40 candidates, re- 
spectively. 

The results obtained in this paper demonstrate the 
usefulness of using multi-band imaging data for a prelimi- 
nary pruning of the cluster candidates and for the selection 
of individual galaxies for follow-up spectroscopic observa- 
tions. As mentioned above we arc currently c;omplcmcnt- 
ing the V- and /-data with /?/? J/sTg-observations which 
will allow us to carry out not only the color analysis as 
presented here but also to estimate photometric redshifts. 
This will significantly improve our ability to confirm clus- 
ters and to select spectroscopic targets. Furthermore, on- 
going spectroscopic observations of EIS clusters selected 
in different redshift ranges will provide key information 
either to confirm the results presented here or to refine 
our selection methods, resolving some of the issues raised 
by the present analysis. 
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